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Abstract: We present a method for the inclusion of finite width effects in the simulation 
of Beyond Standard Model (BSM) physics. In order to test the validity of the method we 
compare our results with matrix elements for a range of production and decay processes 
in the Standard Model, Minimal Supersymmetric Standard Model (MSSM) and Minimal 
Universal Extra Dimensions model (MUED). This procedure has been implemented in the 
Herwig++ event generator and will be available in a forthcoming release. 
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1. Introduction 

It has been believed for some time that there must be new physics at the TeV scale. While 
most theories of physics beyond the Standard Model differ in their internal workings they 
all possess a common feature, the appearance of a spectrum of heavy states that will be 
produced at a collider with sufficient energy. If there is to be any hope of discovering 
these new particles one must have accurate simulations of the signals they will produce, 
and in particular, accurate simulations that can be compared directly with experimental 
data. There are, essentially, two methods of producing events distributed as they would 
be in an experiment. The first is to use a matrix element generator [1-6] to compute the 
full 2 ^ n matrix element for a specific final state of interest, including the decays of 
all of the fundamental particles. This can then be interfaced to a general-purpose event 
generator to produce the particles observed in an experiment. The advantages are that 
the full matrix element calculation gives contributions from non-resonant diagrams which 
may be important and the effects of spin correlations are automatically included. Ideally, 
we would use an n-body matrix element for all processes, however, there are a number of 
issues: 
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1. the time to compute the full matrix element, even if only the resonant diagrams 
are included, grows with the number of final-state particles. This is particularly 
problematic for models such as supersymmetry which can contain long decay chains; 

2. many new physics models introduce a large number of new states, which lead to many 
possible production and decay mechanisms. If the full calculation is used in all cases 
the number of processes required becomes prohibitively large; 

3. any new coloured states will emit QCD radiation in their production and decay which 
must be simulated using the parton shower approximation. 

Given these issues, for the foreseeable future, we will need to use general purpose event 
generators [7-12] that treat the production and decay of heavy particles in a factorised 
approximation. If a generator incorporates the model of interest it can simulate everything 
from the initial hard interaction up to the final colour-singlet particles that interact with 
a detector. The production step is typically a simple 2 — > 2 scattering followed by a series 
of perturbative decays which result in stable particles that are hadronized and decayed to 
give the observed colour-singlet states^. It is vital that these simulations are as accurate 
as possible by, for instance, including effects such as spin correlations [13-15] throughout 
the event simulation. 

Another area where the factorized approach can neglect important physics is in the 
treatment of off-shell effects, as by definition the method must use the narrow width ap- 
proximation to separate production and decay. In the simplest approach all particles are 
on-shell throughout their production and decay [16], which in many scenarios is a good 
enough approximation but in other regimes, e.g. close to thresholds and for resonances, 
such as the boson, where the width can be measured, this approach is unsatisfactory. In 
order to improve the physics of the simulation, off-shell effects must be included in the pro- 
duction and decay stages of the event generation. The aim of this paper is to describe how 
this can be achieved and discuss its implementation in the Herwig-|— |- event generator [7]. 

The next section describes the narrow width approximation in more detail. Sect. ^ 
gives a description of the method by which we include width effects, with specific exam- 
ples for the production and decay stages. Section |^ gives results for various scenarios in 
the Minimal Supersymmetric Standard Model (MSSM) and the Minimal Universal Extra 
Dimensions (MUED) model, some conclusions are drawn in the final section. 

2. Narrow Width Approximation 

In general the evaluation of the matrix element for a process with a high multiplicity final 
state is complex due to the factorial growth in the number of diagrams. The calculation 
can be simplified by using the narrow width approximation, where if: 

1. the resonance has a small width T compared with its pole mass M, T <C M; 

^This is an extremely simplified schematic of event generation. More detail can be found in, for example, 
Ref. [7]. 
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2. we are far from threshold, \fs — M ^ F, where \fs denotes the centre-of-mass energy; 



3. the propagator is separable; 

4. the mass of the parent is much greater than the mass of the decay products; 

5. there are no significant non-resonant contributions; 

the cross section integral can be separated into an on-shell production step followed by a 
series of decays. The separation arises from integrating out the propagators connecting 
each step giving a momentum independent factor 



dg2 



1 



(g2 - M2) + iMV 



(2.1) 



Mr 



If the above assumptions are true one obtains an estimate of the cross section with an error 
of 0(r/M) using this approach. 

In reality, especially when dealing with BSM physics, the approximation is commonly 
used when the above assumptions are not strictly satisfied. While there have recently been 
studies of the validity of narrow width limit for some SUSY scenarios [17-19], nothing 
has been studied in relation to other popular new physics models. In the next section we 
describe how off-shell effects can be simulated with our main focus on BSM studies. In 
Herwig-|— I- the same approach is used for the simulation of off-shell effects in hadron decays, 
allowing the same infrastructure to handle all decays. 

3. Off-Shell Weight Factor 

As discussed above, the narrow width approximation allows the propagator connecting 
production and decay of successive decays to be integrated out. This essentially means that 
part of the phase-space integral is approximated to a constant when the correct assumptions 
are satisfied. Therefore, to improve the accuracy of our simulation we wish to move away 
from the on-shell approximation and include the effects from integrating over the connecting 
propagator. In the past this has been accomplished in a variety of ways. For example, the 
FORTRAN HERWIG [8] program included: 

1. the full three-body matrix element, with an off-shell boson, for top decay; 

2. smearing of fundamental particle masses using a Breit-Wigner distribution; 

3. a more sophisticated Higgs boson lineshape [20]. 

To improve our simulation we include the weight factor (see Sect. ^ for a derivation) 



max 



7r7^2. (m2-M2)2 + m2r2(m)' 



mill 



throughout the production and decay stages, where r(m) is the running width of the 
particle, M is the pole mass and mmin,max are defined such that the maximum deviation 
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Figure 1: Feynman diagrams for tlie process gg ^ q*q wliere tlie Roraan indices give the colour 
representation. 

from the pole mass is a constant times the on-shell width. The weight includes a momentum 
dependence into the calculation of cross sections and decay widths thereby improving the 
approximation to the full matrix element. While for the latter case convoluting the weight 
with the partial width calculation for a particular decay mode is relatively simple, this is 
not the case for the production stage, so we will consider them separately. 

3.1 Off-shell Masses in Particle Production 

For production we need to convolute the weight factor described above with the cross 
section integral. We achieve this by distributing the masses of the outgoing particles 
according to Eq. (3.1). This in itself is a trivial task but what we will show here, with an 



example from supersymmetry, is that gauge invariance can be violated if these masses are 
used naively when calculating the matrix elements. 

Consider the process gg — > q*q, for which the diagrams are shown in Fig. |^, where we 
wish the outgoing squarks to have masses and m4 respectively. Due to the presence of 
external gluons, the Ward identity 

PimP2.-M''"(pi,P2) = (3.2) 

where pi are the momenta of the gluons and M is the total amplitude, must be satisfied. 

After replacing the external polarization vectors with their momenta the amplitudes 
are given by 

Mt = -gl ~ ""^ j [ml - t) t^'d.tt,, (3.3a) 

M^ = {ml - u) 4t, (3.3b) 

A^. = -f (i-n)(4i?c-id.4), (3.3c) 
A^c = f^(44 + 44), (3.3d) 
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Figure 2: The top quark line shape for mP°'° = 174.20 GeV and Ttim^ ) = 1.40 GeV at the ILC 
with the results of Herwig++ for e+e" ^ tt ^ bbW'^W" compared to the Madgraph calculation 
of e+e^ bhW^W" with the Madgraph result including (a) all diagrams containing a top quark 
line and (b) all diagrams excluding those containing a Higgs boson. In both cases the Herwig++ 
result uses our off-shell treatment while the Madgraph result includes all diagrams for the 2^4 
scattering process, including the non-resonant contribution. 



where s, t and u are the Mandelstam variables, are the t- and u-channel masses and 
tfj are the SU(3) colour matrices. Equations (3.3a) and (3.3b) show that for any hope 
of achieving the correct cancellation we must set mt = and rriu = m4. This also 
shows why, even in the on-shell case, one must take care when using widths in scattering 
diagram calculations as these alone can give rise to violations of gauge invariance. The 
total amplitude saturated with the gluon momenta for the g g ^ q*q process is then 



m' 



(3.4) 



so that Eq. (|3.2| ) is only satisfied if = m^. 

This requirement means that, in general, we cannot use off-shell masses when calcu- 
lating matrix elements since if we generate a process such as that shown, we would violate 
gauge invariance. In our procedure the off-shell masses are used when calculating the 
momenta of the outgoing particles involved in the hard interaction but are then rescaled, 
such that ms = m^, for the matrix element calculations. To demonstrate the validity of 
this procedure we compare the line shape of the top quark from Herwig-|— |- and Madgraph 
for the production of a top quark at the ILC, the Tevatron and the LHC. In Herwig-|— |- 
the top quark width is computed using the full three-body matrix element whereas in the 
Madgraph case just the two-body decay of the top quark is used, due to rapid growth in 
the number of diagrams that are required. In all cases Herwig-|— |- generates the 2 — > 2 
production process for the tt pair followed by the three-body decay of the top quark using 
the treatment of off-shell effects described in the text. Madgraph was used to calculate 
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Figure 3: The top quark line shape for m^"^" = 174.2 GeV and Ttim^"^") = 1.4 GeV at (a) the 
Tevatron and (b) the LHC. The black line denotes the results from Herwig++ with the outgoing 
masses rescaled to their on-shell value, the dot-dash line denotes Herwig-|-+ with the outgoing 
masses rescaled to their average value and the crosses denote the Madgraph results. The Herwig++ 
results were generated using the 2^2 production process for ti followed by the three-body decay of 
the top quarks. The Madgraph results use the full matrix clement for the production of bbW^W^ 
to order OgO^, including all diagrams, both resonant and non-rcsonant diagrams, containing a top 
quark line. 



the 2 — > 4 matrix element for the production of 6 6 including the non-resonant 

diagrams. 

To ensure that the amplitudes generated by Madgraph [1] were gauge invariant, the 
"fudge-factor" scheme [21] was used. This involves calculating the full amplitude with- 
out the inclusion of the width for any off-shell propagators and then multiplying the full 
amplitude, including non-resonant contributions, by 

p2_M'i + iMT 

for any propagator that can be on-shell, i.e. for which it is possible for = within the 
physically allowed phase space. This is the simplest approach that ensures the amplitude 
is gauge invariant [21], although it has the drawback that the non-resonant diagrams are 
affected. A more detailed discussion of the issue of gauge-invariance when including non- 
resonant diagrams can be found in [21]. 

For the ILC case, Fig. |2|, the Madgraph result is shown for both the process including 
only the diagrams with a top quark line and also the process including all electroweak 
diagrams, resonant and non-resonant, excluding the Higgs. There is excellent agreement 
between our results and those performed with the full matrix element giving us confidence 
in our procedure. For the hadron colliders we must consider the rescaling since there will be 
processes such as gg ti that will, as discussed above, violate gauge invariance when we 
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Figure 4: The line shapes for (a) the SUSY partner to the u-type doublet field and (b) the level-1 
KK partner to the u-type doublet field. 



take the top quark off-shell. Here we will compare two choices for the momenta rescaling, 
first rescaling such that the masses have their on-shell value and second rescaling to the 
average value of the outgoing masses (ms + m4)/2. The results for the Tevatron and the 
LHC are shown in Fig. ^. The Tevatron results are in excellent agreement with the matrix 
element for both choices of rescaling and the LHC is good agreement except for the tail 
where there is a small deviation. It is clear that either choice for the value of the rescaled 
mass gives good agreement with the matrix element but in the LHC case choosing to rescale 
to the average value of the outgoing masses gives slightly better agreement with the full 
calculation. 

Figure ^ shows the mass distributions for a left-handed up squark in the MSSM^ 
and the KK-partner of the doublet quark in the MUED model respectively. The mass 
spectrum for the MUED case is matched to the SUSY spectrum at SPS point 2 [22] where 
rriu^ = 1560.97 GeV, Tmssm = 70.22 GeV and Tmued = 312.76 GeV. 

This example is at the extreme of where this method should be applied since, especially 
in the MUED case, the width is large and in general there could be sizeable non-resonant 
contributions. 

3.2 OfT-Shell Effects in Particle Decay 

Many new physics models have spectra that result in long chains between the production 
of a resonance and a stable state. As mentioned previously the simplest approach in 
dealing with these chains is via a series of on-shell cascade decays. While this may be an 
appropriate approximation in some kinematic regions, in others, i.e. when the decaying 

^It is technically incorrect to say that a scalar particle has a helicity state but the terminology allows 
for easier distinction between the partners of the SU(2) doublet and singlet quarks. 
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particle is close to threshold, the effects from the off-shell propagator must be taken into 
account. 

This can be achieved by including the weight factor from Eq. (|3.1| ) in the calcula- 
tion of the partial width of a selected decay mode. For example, consider the decay 
g — > 6 61 ^ X2 ^ ^! the partial width is 



^{~g^xlbh) = - dm^ 1 r(g^66i , (3.6) 

where the widths inside the integral are evaluated for the off-shell mass m. The limits on 
the integration are determined by the on-shell width and are set such that the maximum 
deviation from the pole mass of 61 is ST. As the intermediate particle is a scalar, the inclu- 
sion of the weight factor should give exact agreement with the full three-body calculation 
providing the integral is performed over the same phase space. Fig ^ demonstrates this 
for SPS point la where the three-body phase-space is restricted to the same as the two- 
body case. The spectrum was produced using SPheno 2.2.3 [23] where m-^^ = 515.27 GeV, 
F(m^J = 3.83 GeV and m^o = 180.58 GeV. The mass of the 6-quark is mi, = 4.20 GeV, 
which is the default value in Herwig++. The on-shell result is also included for reference. 

The agreement between the full matrix element calculation and our results show that 
the approximation is valid. 



4. Examples 

Here we present a range of processes in the MSSM and the MUED model demonstrating 
the consistency of the inclusion of off-shell effects in Herwig-|-+. In the decay examples the 
comparison is always to the full three-body result, which was included in Herwig++ as a 
modification to the current released version and will be included in a future version along 
with the simulation of finite width effects. 
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4.1 Decay via an Off-Shell Fermion 

A possible two-body decay of the ii squark in the MSSM is ii — > Xi^- If i^ti ~ "^x? ~'~ 
then the effect of the width of the top quark must be considered. We choose the decay mode 
ii ^Xit^Xi W+b at SPS point la [22] where m-o = 97.04 GeV with mj = 174.20 GeV, 
mw = 80.40 GeV and nn, = 4.20 GeV. The threshold values for the on-shell two- and 
three-body decays of the ii are 271.24 GeV and 181.64 GeV respectively. Figure ^ shows 
partial width of the ti as a function of its mass for the three-body, two-body off-shell and 
two-body on-shell results. 

Unlike the gluino decay example in Sect. 3.2 the boson, like the top quark, has a 



measured decay width and this should be treated properly. In the example shown in Fig. ^ 
the running width for the top quark is calculated from its full three-body matrix element 
to a 6-quark and a pair of light fermions which includes the full effects of the width. 
The agreement between the two-body off-shell and three-body results shows that this is a 
valid approximation to use. Also, despite the extra factor of {ji+m) in the numerator of 
the fermion propagator, there is still good agreement between the full three-body result 
and the two-body result with weight factor even though the factor does not attempt to 
include this. This is due to the numerator factor being largely responsible for propagating 
spin information rather than altering the kinematics. 

4.2 Decay via an Off-Shell Gauge Boson 

In the MSSM there is a coupling between the boson and the gaugino sector allowing for a 
three-body decay of the second neutralino to a pair of light fermions via an intermediate 
boson. The presence of a spin-1 rather than a spin-0 particle alters the form of the partial 
width as the decay is now j?-wave and not s-wave, as in Fig. |. To illustrate that the weight 
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Figure 7: The partial width for X2 ^ Xi -^^ ~* bbxi in the MSSM. The lower panel gives the 
value of (rthrco/roff - !)■ 

formula works just as well in this situation we choose the decay chain X2 ^ Xi^^ ^ Xibb 
at SPS point la where m^o = 97.04 GeV with Mz = 91.19 GeV and mb = 4.20 GeV. 

Figure ^ shows the results for the above decay and demonstrates that there is good 
agreement between the full three-body result and the two-body approximation for an in- 
termediate vector particle. Another example of a possible p-wave decay is u* ue\~ e'^ in 
the MUED model where the intermediate particle is the level-1 KK-Z" boson. For param- 
eter values = 500 GeV and h.R = 20 the relevant masses are M^o = 535.81 GeV and 

= 504.25 GeV. The partial width is shown in Fig. ^, again with both the three-body 
result and two-body via an off-shell Z^. 

5. Off-Shell Cross Sections 

In the narrow width approximation, the cross section for a particular final state is com- 
puted by taking the on-shell production cross section to an intermediate resonance and 
multiplying by the branching fraction to the final state of interest. If the on-shell mass of 
the resonance is close the threshold for the decay into the final state then the narrow width 
approximation is invalid and one should calculate the full matrix element. As described 
in Sect. ^ and demonstrated in Sect. |^ we can include a weight factor in production and 
decay to simulate off-shell behaviour. In the case of calculating cross sections for specific 
processes, this amounts to including the effect of propagator widths in the Monte Carlo 
estimate of such a quantity. It is important to note, however, that a general purpose event 
generator that starts from a 2 ^ 2 hard scattering and then perturbatively decays the pro- 
duced resonances can never include non-resonant contributions. This is a fundamental limit 
of the approximations used to generate the events. Nevertheless, a good approximation 
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Figure 8: The partial width for u* — > u e° e+ in the MUED model. The lower panel gives the 
value of (rthrcc/roff - 1). 



can still be achieved providing one uses the simulations with care^. 

An example of a process that has no non-resonant contributions is the production of 
a strange squark via uJ — > g xt sls, the diagrams for which are shown in Fig. |9|. 
The results for the ratio of the off-shell to the on-shell cross section as a function of the 



strange squark mass are shown in Fig. IC for SPS point la. The ratio is constant, with the 
off-shell result smaller due to the integration limits no longer being taken to infinity, until 
rrigj^ « O.Srrig where we are in the threshold region for the decay of the gluino. The sudden 
steep rise as the mass ratio approaches unity is due to the on-shell cross section going to 
zero at threshold. 

There is a counterpart process to that in Fig. ^ for MUED, where the xt replaced by 
the boson, the g by the gi and sl by the s*. The ratio of the on-shell to the off-shell 



cross section for this process is also shown in Fig. 10 where the masses for the MUED 
particles have been matched to SUSY spectrum to give a fair comparison. It is apparent 
here that the spins of the underlying model play only a small role in determining the value 
of this ratio as the results are similar and while the absolute values of the cross sections 
may differ greatly, taking the colour octet object off-shell affects only the kinematics. 



6. Merging Two- and Three-Body Decays 

In Sect. we demonstrated the accuracy of including an off-shell weight factor by compar- 
ison with the full three-body matrix element for a variety of processes. For each process 
considered the width was plotted over the entire kinematic range, rather than restricting 
to the region where the decay would be applicable, to give a full comparison. In a real 



^In some specific cases non-resonant effects can be modelled by using a modified form of Eq. (3.1), for 
example the Higgs lineshape [20]. 
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Figure 9: The diagrams contributing to the process u d - 
to a strange quark and a left-handed strange squark. 
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Figure 10: The ratio of the off-shell and on-shell cross section for the process ud ^ xt ^ls (black 
line) and its MUED counterpart (dashed line). 



simulation there is a choice over which point we should change between using a two- and 
three-body decay of the particle. If both decays were treated on-shell then the point would 
simply be the threshold of the two-body decay but when including of off-shell effects for 
the two-body decay, the choice is not so simple. 

Here we use the three-body decay X2 ~^ Xi e"*" e~ in the MSSM to study this effect. 
The full three-body decay is mediated by a boson and both the left- and right-handed 
selectron, giving an interference between the different channels. If the decay occurs as a 
series of cascades with a weight factor then these interference effects will be neglected. To 
judge the extent of the interference we compare the partial width for the decay via the 
full three-body matrix element and the three-body matrix element with the diagram 
removed and performed cascade decay. 



Figure 11 shows the results for a range of selectron masses, where both cl and cr 
are degenerate, with M^o = 120.00 GeV and Mz = 91.19 GeV. For a sufficiently large 
selectron mass there is good agreement between the two methods as there is only a small 
interference with the boson diagram. However, as the mass is lowered so that the decay 
of the X2 through the selectron mode becomes closer to being on-shell, the interference 
effects, particularly just above the Z^ threshold, become significant and the full three- 
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Figure 11: The partial width for the decay X2 ~^ Xi e~ where the selectron masses are indicated 
on the plot. The solid line is the full three-body partial width and the crosses are for the three-body 
decay with the diagram removed plus r(x2 Xi Z^) ^ BR{Z'^ — > e+ e^). The dashed line shows 
the two-body on-shell cascade. 



body calculation is necessary in this region. For the final case, Fig. 11c, where the selectron 
modes are on-shell for the whole range there is quite different behaviour. The partial width 
now smoothly passes over the threshold and there is exact agreement with the full three- 
body result indicating that there is very little or no interference. Given these results it 
seems reasonable to use the threshold of the on-shell two-body decay as the point where 
the change from a three-body to a two-body decay with weight factor occurs. 



7. Summary 

Given that we do not know what type of new physics will be discovered at future colliders, 
it is necessary to have access to accurate simulations of these new models. An important 
consideration when studying new physics scenarios within a Monte Carlo event generator 
is the simulation of off-shell effects. Here we have demonstrated a consistent algorithm 
for their inclusion in a general-purpose event generator using a variety of processes from 
the Minimal Supersymmetric Standard and Minimal Universal Extra Dimension Models. 
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Figure 12: A three-body decay consisting entirely of scalar particles where the external particles 
a, b, d and e are all on-shell. 

The approach described here has been implemented in the Herwig-|— |- program and will be 
incorporated in to a future release. 
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A. Derivation of Weight Factor 

The weight factor introduced in Eq. (|3.lD can be derived by considering a three-body decay 



that consists entirely of scalar particles. Using the notation in Fig. 12, the decay rate is 
given by 

T{a^b,d,e) = J d(f>3 {pa;Pb,Pd,Pe)^\M3,\ (A.l) 



where dff)^ is the three-body phase-space and jA^sl is the spin-averaged matrix element. 
The phase space can be written recursively as [24, 25] 

d(p3{pa;Pb,Pd,Pe) = d(p2{pa;Pb,q)d(p2{q;pd,Pe){'2T^fdq^, (A.2) 

where d<j)2 is a two-body phase-space factor and q is the momentum of the intermediate. 

For the matrix element assume that the intermediate particle of mass M has n two- 
body decay modes so that 



-2 



I A yl I 



where \A4a\ and jA^cJ are the two-body spin-averaged matrix elements and T{q) is the 
width of the intermediate at scale q. Substituting Eqns. ( A.2 , A.3| ) into Eq. ( [A.l| ) gives 



ma J {q^ - M^) + q^T^{q) J J ^1 ^ 

(A.4) 
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The third integral in Eq. ( |A.4| ) can be recognised as r(c — > d, e) using 



r(o b, c) 



2mn 



\Ma\ d02 



and the second as T(a b,c), with the intermediate at scale q, giving 

(g2 -M2)2 + g2r2(g)' 



T{a^b,d,e) = ^ j dq^Ta 



(A.5) 



(A.6) 



The weight factor is then identified as 



w = — f dq^ g 



(A.7) 



The case for production followed by decay follows similar arguments and the same factor 
is found. 

It should be noted from the discussion in Ref. [20] that use of the running width in 
the propagator of Eqn. ( [A.4| ) is only valid if T{q) ~ q for large q. If T{q) were to grow 
faster than this then the extra terms are dominant and the propagator becomes of 0(1 /a), 
the coupling of the c ^ d,e decay. If T{q) grows linearly with q the extra terms are just 
unenhanced higher order corrections. 
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